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Abstract 

We begin with a prior observation by one of us that Thomas preces- 
sion in the nonrelativistic hmit of the Dirac equation may be attributed 
to a nonabehan Berry vector potential. We ask what object produces 
the nonabehan potential in parameter space, in the same sense that the 
abelian vector potential arising in the adiabatic transport of a nonde- 
generate level is produced by a monopole, (centered at the point where 
the level becomes degenerate with another), as shown by Berry. We 
find that it is a meron, an object in four euclidean dimensions with 
instanton number i, centered at the point where the doubly degener- 
ate positive and negative energy levels of the Dirac equation become 
fourfold degenerate. 
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Consider the adiabatic transport of a nondegenerate level of a hermitian hamilto- 
nian around a closed loop in parameter space. If \x{R)) is the instantaneous eigenket 
of the hamiltonian at a point R in parameter space, Berry showed that 

(i) the state vector picks up a phase factor due to a vector potential 

A, = {x\d,x) (1) 

in addition to the integral of the energy along the circuit and that 

(ii) The source of this potential is a quantized magnetic monopole located at a point 
where the level becomes degenerate with another. In other words, the field tensor 
Ff^i, derived from is that of a monopole. 

The proper mathematical framework for describing Berry's work was established 
by Simon. 

It was then pointed out by Wilczek and Zee that if the level being transported 
has an m-fold degeneracy, the effect of a closed circuit will now be given by U{m) 
matrix obtained by a path-ordered integral of a nonabelian vector potential. These 
authors provided many illustrations. 

The references mentioned up to this point, along with many others and very useful 
commentary may be found in the collection edited by Shapere and Wilczek 

The problem of two-fold degenerate levels becoming four-fold degenerate was dis- 
cussed in its generality by Avron ei a/ , using spin systems as a model. The relation 
of our work to theirs will be discussed at the end. 

Our problem here stems from the observation by one of us that Thomas pre- 
cession of electron spin in the nonrelativistic limit could be understood as arising due 
to nonabelian SU(2) Berry potential. The idea is as follows. For the free Dirac equa- 
tion, the four levels at each momentum break up into two (spin) degenerate pairs with 
equal and opposite energies. Since there is a large gap separating the positive and 
negative levels, they will not mix under adiabatic evolution. (In practice the coulomb 
potential of the nucleus causes the adiabatic evolution). It was then a straightforward 
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matter to compute the vector potential 



A 



ab 



(a I (9^6) 



(2) 



starting from the degenerate positive energy spinors labeled by a and b and to show 
that the effect of this potential was to produce the Thomas precession. (It is of course 
possible to do the same for the negative energy states. Throughout this paper we will 
focus on just the positive energy states. ) 

Here we ask and answer the following question: What is the source of this non- 
abelian potential? We find it is meron, an object with Pontrayagin index equal to |. 
It is not necessary to be familiar with merons to follow this paper since their relevant 
features will be described later on. 

Let us begin with the Dirac hamiltonian 



where the three 7's and (3 are 4x4 matrices which anticommute and whose square 
is unity. Introduce 4 euclidean gamma matrices 



H — 'y ■ p + (3m 



(3) 



7/. = (7,/3)) 



(4) 



obeying 



(5) 



and a four vector 



(6) 



Note that the fourth component of p^ is m and not the energy. Then 



H^^ = ^^p^ (// = 1,2,3,0). 



(7) 



We shall use the representation 
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7m 



< 



(8) 



where 



= (<7, il) 

and cr are the Pauh matrices and / is the 2x2 identity. 
Since 

p^p^ = p2 ^ p2 + ^2 



(9) 



(10) 



it follows the eigenvalues of H are — ±p. Since there should be a total 

of four levels, it follows they come in doubly degenerate pairs. The degeneracy may 
be traced back to fact that 



H = 727o-f^*7o72 



(11) 



which is Kramers' degeneracy. The fact that the levels come with equal and opposite 
energies is due to the fact that 



75 = 7i 727370 = 

anticommutes with H and is unitary. 
Let us now introduce a unitary matrix 



/ 
-I 



(12) 



(T ■ p — im 

P 



(13) 



in terms of which 

pW 
pU 

and the positive energy eigenvectors are 



H 



(14) 
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UXa 



(15) 



where Xa and Xb are any two orthonormal two-component spinors. We choose them 
to be the canonical basis 



The nonabehan vector potential is now given by a 2 x 2 matrix 



Since x is a canonical basis vector we can write in compact notation 



Note that due to the factor of 1/2, this is not a pure gauge. 
Explicit computation shows that 



Ao = ^^-p A 



— {ma + (TX p) 



The vector potential has no radial component: 



The next step is to compute the field strength 



Given eg. (PI^D , it follows that 



(16) 



(17) 



(19) 



(20) 



(21) 



(22) 



Although this expression may be evaluated in a straightforward but tedious man- 
ner, it is instructive to employ some tricks analogous to the ones used by Berry in 
the abelian case. 
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Let us first introduce some conventions. The kets \a), \b), \c) will stand for positive 
energy eigenvectors while 1^4), |C) will stand for negative energy eigenvectors. We 
need two lemmas: 
Lemma 1: 

(MB) ^ m (23) 

proven by taking the derivative of {a\H\B) — and using {da\B) = —{a\dB). 
Lemma 2: 

{a\dH\b) = a^h (24) 

proven by taking the derivative of {a\H\b) = for the case a ^ b. 
Now we have 

F^t^d^dub)- {fi^u) + {a\d^c){c\d,b)- {fx ^ u) (25) 
^ {d^a\d^b) - {ix^ v) - {d^^a\c){c\d^h) - {n ^ v) (26) 
^ {d^a\A){A\d,\b) - ifx^u) (27) 

where in the last line we note that the intermediate states are restricted to those of 
negative energy. We now invoke Lemma 1 to write 

F^',^^{a\d^H\A){A\d,H\b)- ix^v (28) 

We now argue that we can enlarge the set of intermediate states to include the positive 
energy states well. First let a ^ b. Then no matter which state c we introduce, it 
will be unequal to either a or b, and the term will vanish by Lemma 2. Finally if 
a — b, then the contribution from c will vanish if c a thanks to Lemma 2 and by 
antisymmetry if a = 6 = c. 

We can now use completeness to reach the nice result 

Fi-^2i^\i^H^d'^H]\b) (29) 
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(30) 
(31) 



where it is understood [i^^ v and where 



a^^ai. 



(32) 



If /i and p are spatial indices, 7^jy is self-dual (same in the upper and lower blocks 
which are eigenspaces of 75 with eigenvalue ±1), and anti-self-dual (of opposite sign 
in the two blocks) if one of the indices is 0. 

Note however that the F^^, for our problem is not self-dual or anti-self-dual since 
it is given by the projection of 7^1, into the positive energy subspace. Indeed since 75 
anticommutes with if, its eigenvectors cannot be also eigenvectors of H . 

Using the explicit formulae for the eigenvectors, given by eqs,(|l5]-0) we find that 

= 4^ [^'^^'^ + U^^l^A ■ (33) 

In the problem of Thomas precession, the particle goes around in a tiny loop (in 
the nonrelativistic limit) in some plane, say the the x — y plane, at fixed m. The 
relevant field is the "magnetic field " along the z-axis, given in this limit by 



F = B 



2ni? 



(34) 



in agreement with Reference 6. (In that paper the vector potential was given in a 
different gauge due to a different choice of eigenvectors.) 

But now that we had the field in all of parameter space (and not just the nonrel- 
ativistic region) and we decided to explore some of its properties. We found that 



D^F^u = d^F^, + [A^, F^,] = DF = 
Tr{FF) = -8tt''6^^\p). 



(35) 
(36) 
(37) 
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The first equations has a counterpart in Berry's abehan monopole: there the 
magnetic field is radial, R x B = 0, which can be written as RiFij = 0. It is possible 
to choose the vector potential so that it too has no radial component. Likewise in 
our problem, only the circulation in tangent planes is nonzero. 

The second equation tells us that F solves the Euclidean Yang-Mills equations. 
It is however not self-dual or anti-self- dual, nor is it required to be, since the corre- 
sponding action (the integral of TrF"^) is logarithmically divergent. 

The last equation, ( which refiects the fact that Tr E ■ B = here) tells us that 
the instanton density is zero everywhere, except possibly at the origin where the 
field strength is singular. Indeed one can argue that there must be a delta-function 
singularity there as follows. 

The instanton density Q may written as the divergence of a vector as follows: 
(See for example Reference 7 for details and references) 

Q = -^Jr{FF) (38) 

= d^K^ (39) 

1 2 

= -^^i.vapTr[A,{d^Ap + -A^Ap)] (40) 

Thus the total instanton number inside any volume can be found by doing the surface 
integral of K^. Since must have the form cp^/p'^, by dimensional analysis, we can 
find c by considering a point on the m-axis. It is readily seen that c = l/47r^ which 
fixes 

Tr{FF) = -M^5'^^\p) (41) 



which in turn implies Q = 1/2 as per eq.( 

At this point it became clear that the field was that of a meron. The meron 
was first discovered as a solution to Yang-Mills equations by De Alfaro et al 
with instanton number 1/2 and invoked by Callan et a/ Q as a configuration that 
could enter the functional integral for the Yang-Mills field and explain confinement by 
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producing the area law decay for the Wilson loop. We paraphrase their description of 
the relation between the meron and the instanton, which we found very enlightening. 
We will however do it in terms of the two-dimensional instanton which arises in the 
0(3) sigma model since it is easier to visuahze. Consider the instanton of unit size 
which is obtained by placing a sphere of unit diameter on top of the plane and doing 
a stereographic projection and assigning to each point on the plane the unit vector 
associated with its inverse image: thus the origin gets assigned the south pole, the unit 
vector slowly starts rotating upwards as we move out and the equator gets mapped 
into a unit circle . As we move along this circle, the field lies in the plane and is 
radial. As we move further out, the field tilts up even more and finally at infinity, 
points upwards. If we measure the instanton number enclosed within any circle, it 
will grow from zero to a half by the time we reach the unit circle, and reach unity 
by the time we integrate out to infinity. Now we deform the map as follows. First 
we scale down the region inside the unit circle, which carries the lower hemisphere, 
to a circle of radius r < 1; push the region outside the unit circle (which carries 
the upper hemisphere) to outside a radius R > 1 and assign to the entire annulus 
r < 1 < R, the radial equatorial field. The instanton density reaches the value 1/2 
at r, stays fixed out to it! (since the entire annulus maps back to just the equator) 
and starts growing to unity as we go past R to infinity. // we now let r ^ and 
R ^ oo we get the meron. It will look just like the vortex of the x-y model, have 
a logarithmically large action, (In i?/r) and carry no instanton density except at the 
origin. Our meron is related to the four dimensional instanton in just the same way: 
one takes the instanton, squeezes the region with half the winding number into the 
origin, sends the other half out to infinity and fills all of space with a configuration 
that has zero Tr FF. 

We now relate our work to the general analysis of Avron et al of hamiltonians 
with Kramers degeneracy, brought to our notice as this work was nearing completion. 
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These authors showed that the generic hamiltonian with Kramers degeneracy hves 
in a five dimensional Chfford algebra. Let us call the coefficients = (p, m, M). 
This corresponds to our adding a term M75 to our H. Because of our interest in 
the physical problem of the Thomas precession of the Dirac electron, we shced this 
space along the plane M = 0. The relation of our four dimensional analysis to the 
five dimensional one is best explained by a simpler analogy. Consider the monopole 
that arises in Berry's treatment of the abelian problem. It lives in three dimensions. 
It can be described by a vector potential with no radial component and a field tensor 
which too has circulation only in the planes tangential to the position vector. By 
surrounding the monopole with a unit sphere S'^ and measuring the flux (in proper 
units) we can obtain a topological invariant of the map (first Chern number) from 
the space of states to S"^. On this sphere the field strength is uniform and finite. 
Suppose one were now to slice this space along the plane z = 0. On this plane one 
will see a radial field of the x-y vortex, i.e., the meron. The unit sphere S"^ will be 
sliced along the unit circle enclosing the meron. Consider now the integral of the 
monopole vector potential = (xl^^^x) along this circle. By Stokes' theorem this 
equals the fiux enclosed in any surface bounded by it. If one uses as the surface the 
upper hemisphere of the 5"^, one gets half the topological charge, the contribution 
being uniform on the hemisphere. The other choice for the surface is the interior of 
the circle in the x-y plane, i.e., the disc with as its boundary. In the latter case the 
fiux is zero except at the origin where the contribution is 1/2 due to all the fiux fiowing 
upwards into the upper hemisphere. Thus the person restricting himself to the plane 
can reconstruct the topological index associated with the S"^ by doing an integral on 
the he sees, thanks to Stokes theorem. (He must then remember to double to the 
answer for the lower hemisphere in this symmetric problem, or more generally, take 
the difference with another oppositely oriented line integral of the vector potential 
that describes the lower hemisphere.) 
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Returning to our problem, we must begin with H = ]9^7^ = in five dimen- 
sions, find its eigenvectors and define = {x\d^x) ^-^id find the ten component field 
strength JF^,^. The result will be just as in eq. (^1]) (thanks to the lemmas) 

= ^((^Mb) (42) 

with the obvious change in the definition of p and the range of indices. The field is 
still singular at the origin in five dimensions, obeys p^J^^v = and DjF = 0. The 
former implies that the field is tangential. (This is true in any gauge. Also one can 
pick a gauge, as we did, in which itself has no radial component.) Thus at each 
point Pfj^, T has only six nonzero components with "transverse" indices. We shall 
refer to them as F . This configuration is thus a natural generalization of the Berry 
monopole configuration in three dimensions. 

Avron et al show that if we surround the origin by a unit sphere S^, 

• The (suitably normalized) integral of Tr FF on the sphere (where the dual is 
with respect to the epsilon symbol Pxexf^uafs) is ±1 and measures a topological 
invariant (second Chern number or Pontrayagin index) of the map from space 
of states to S'^. 

• F is is self -dual or anti-self-dual depending on the sign of the energy. 

Since all points on the sphere are equal by symmetry, let us go to one with just M = 1, 
rest equal to zero. Now p\e\^yai3 becomes just the six nonzero components of 

JF have their indices going from to 3. Since H = fi'j^ here, the eigenstates of H are 
indeed eigenstates of 75 and the field is self(anti)-self-dual as observed earlier with 
respect to Sf^^a/s- 

By focusing on the Thomas precession of the physical Dirac equation, we sliced 
this space along M = 0. Our plane contains the singularity at the origin and the 
sphere which is the equator of the S^. The integral on this sphere of i^T^, defined 
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in eg ■ (^OD , gives, by Stokes theorem, either half the topological charge associated with 
the S*^ (if we view the 5''^ as the boundary of the hemisphere of S^) or the meron 
charge if we view it as the boundary of the disc (or ball) contained in the slice. (In 
the language of forms, since Tr FF is closed on the 5*^, it can be can be written 
locally as the derivative of a three form, which is just the dual of K^j^. The integral 
of the topological density associated with the upper hemisphere can be written via 
Stokes' theorem as the integral over the equatorial S*^ of this three form. But by 
Stokes theorem, this also equals the integral of Tr FF within the ball in the M = 
slice bounded by the same S*^, which gives now the meron's winding number of a 
half.) 

To summarize, we looked at the adiabatic evolution of the Dirac hamiltonian and 
asked what singularity in parameter space produces the nonabelian vector potential 
which is behind the Thomas precession and found that it was the meron. Although we 
computed the connection for all possible transports in parameter space, only closed 
paths in momentum at fixed m arise in atomic physics since m is invariant. On the 
other hand, in cosmological models in which m can vary (say with some Higgs field, 
which gives the fermion its mass) other components of the field can produce observable 
effects. The connection to merons is amusing and possibly has other implications. 
Callan et al offer many arguments for the meron's previleged role before proposing it 
as a configuration that can make important contributions to the Yang-Mills functional 
integral. That the meron appears naturally in the Berry analysis provides yet another 
argument. As for the topological ideas discussed here, they are naturally subsumed 
by the very general five dimensional analysis of Avron et al . However the four 
dimensional slice produces singular configurations interesting in their own right, just 
like the x-y model vortices that arise on slicing the monople. We hope readers will find 
the concrete example of the Dirac problem (with its relation to Thomas precession) 
and the relation between the S"^ instanton of Avron et al in five dimensions and the 
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meron that appears in our four dimensional slice, a useful addition to our knowledge 
of instantons, merons, Thomas precession and the Berry connection. 

RS acknowledges some wonderful discussions with Greg Moore. This report was 
supported by an NSF Grant DMR 9120525. 
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